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Abstract 

A comparison of specific interplanetary conditions for 798 magnetic storms with 
Dst < -50 nT for the period 1976-2000 was made on the basis of the OMNI archive 
data. We categorized various large-scale types of solar wind as interplanetary drivers 
of storms: corotating interaction region (CIR), Sheath, interplanetary CME (ICME) 
including magnetic cloud (MC) and Ejecta, separately MC and Ejecta, and "Inde- 
terminate" type. The data processing was carried out by the method of double 
superposed epoch analysis which uses two reference times (onset of storm and the 
minimum Dst index) and makes a re-scaling of main phase of storm in a such way 
that all storms have equal durations of the main phase in the new time reference 
frame. This method reproduced some well-known results and allowed us to ob- 
tain some new results. Specifically, obtained results demonstrate high importance 
of Sheath in generation of magnetic storms as well as a significant differences in 
properties of MC and Ejecta and in their geoeffectiveness. 
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1 Introduction 



One of the important aims of Space Weather is investigation of interplanetary 
drivers of magnetic storms. It has been well known for a long time that the 
most important parameter leading to geomagnetic disturbances and, in par- 
tially, to magnetic storm generation is negative (southward) Bz component of 

interplanetary magnetic fi e ld (IMF! (IFairfield and Cahilll . 1966 ; iRostoker and Falthammaii . 
19671 : iRussell et all Il974l : iBurton et all Il975l : lAkasoful . Il98l[ ). Because IMF 
lies in ecliptic plane under steady interplanetary conditions and substantial 
52; < is observed only in disturbed types of solar wind (SW), it was found 
in many investigations that interplanetary coronal mass ejections (ICME) and 
corotating interaction regions (CIR) are the most important drivers of mag- 
netic disturbances on the Earth. Therefore, it is natural to categorize these so- 
lar wind drivers during a stu dy of magnetic storm generation (see reviews and 
recent papers, for instance. (ITsurutani and Gonzalezl. Il997 



1999 



2006 



Gonzalez et al 



Yermolaev and Yermolae v. 2002: lHuttunen and Koskinen 



2004 



Borovskv and DentonU2006: Denton et a l. I I2OO6I; iHuttunen et al.l 



Yermolaev et al.l . 



2006 



Yermolaev et al., . ,2007a .b.c: Pulkkinen et al.l . l2007a .b) and references therein) . 
Nevertheless, the question what interplanetary conditions results in specific 
magnetic storm features is still open and it does not allow us to precisely pre- 
dict reaction of the magnetosphere. The progress in solving this problem may 
be connected with development of methodical approaches. 



Usage of the peak-to-peak method when minimum of IMF Bz < (or elec- 
tric field Ey = Vx * Bz) was compared with extreme values of the Dst and 
Kp indexes did not allow one to find significant differences between these in- 
terplanetary and magnetospheric param eters for magneti c storms generated 



by different drivers (see, for instance, (lYermolaev et al.l . l2007d ) and refer- 



ences therein). Analysis of time evolution in interplanetary parameters us- 
ing the superposed epoch analysis (SEA) method was more informative and 
showed some characteristic conditions during storms and several differences 
in interplanetary parameters for different drivers (see Table 1). Although 
the CIR, ICME and Sheath (compressed region before ICME) as important 
driver s of magnetic storm have been discussed in the literature for a long 

time (Gosling et al. ^ 1991 ;^ Tsurutani and Gonzalezl. 19971; Vieira et al. . 2004 ; 

Huttunen and Koskinen . 20041: Yermolaev et al. . I2OO5I : lYermolaev and Yermolaevl. 



20061 : lAlves et al.l . l2006l : iBorovsky and Dentonl . |2006|) , still some authors did 



not separate large-scale types of solar wind for superposed epoch analysis 



(iDenton et al.l . l2005l : IZhang et al..l2006l:lLiemohn et al. 



20081: llhe et al.l.l2008h 



2006 


). Only recent papers analyzed separately CIR, Sheath and body of ICME 


(Huttunen and Koskinen, 12004: 


Yermolaev et al.. 


2006; 


Huttunen et al., 2006; 


Yermolaev et al.. 


2007a 


be: 


Pulkkinen et al.. 


2007a). The choice of zero (ref- 



erence) time for SEA is important and substantially influences on results 
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( lYermolaev et al.l . |2006| . l2007al : lllie et al.l . |2008| ). In the most part of previ- 
ous papers the authors used the peak of Dst index as a zero t ime for SEA 
( jMaltsev et al.l . 1 19961 : iLoewe and Prolsd . 119971 : IZhang et al.l . |2006|) . This choice 
is convenient for studying the end of main phase and beginning of recov- 
ery ph ase of storms because duration of the raain phase lasts from 2 to 15 



hours (I Vi chare et al.l . l2005l : iGonzalez and Echerl . l2005l : lYermolaev et al.l . 12006 



2007bl ). inside of an interval with duration of several hours the parameters 
measured before and after onset are averaged simultaneously and specific con- 
ditions resulting in storm onset cannot be studied. The onset time as zero 
time of S EA allows one to i nvest i gate in terp lanetary sources and in itial part 



of storms (lYermolaev et al.l . |2006| . l2007al Jb|]d: iPulkkinen et al.l . l2007al ) . For ex- 



ample, this analysis showed that storms initiated by Sheath have sharper and 
shorter main phase than storms initiated by an other interplanetary drivers 
(lYermolaev et al.l . l2007bl : IPulkkinen et al.l . l2007ah . 



In this paper we study interplanetary conditions resulting in magnetic storms 
on the basis of the OMNI database during 1976-2000 and categorize 6 types of 
solar wind: (1) CIR, (2) ICME, (3,4) two types of ICME - Magnetic cloud (MC) 
and Ejecta, (5) Sheath before MC and Ejecta, and (6) "Indeterminate" type 
(OMNI does not contain sufficient information for identification of the type). 
We use "double" (with two reference times) ESA method, that is, we re-scale 
duration of main phase of all storms in a such manner that, respectively, onsets 
and minima of Dst index for all storms coincide and study interplanetary 
conditions leading to start and end of main phase of magnetic storms induced 
by these 6 interplanetary drivers. 



2 Method 



The basis of our inve stigation is 1-h interpl a netar y and magnet ospheric data 



of OMNI database (iKing and Papitashvili |200J). We made our own data 
archive including OMNI data and calculated (using OMNI data) additional 
parameters. Using level criteria for key parameters of SW and IMF (veloc- 
ity, temperature, density, ratio of thermal to magnetic pressure, magnitude 
and orientation of magnetic field etc.) we defined corresponding large-scale 
types of SW and possible err or of this identification for every 1-h point of 
the archive during 1976-2000 (lYermolaev et al.l . l2009al ). Our identification of 
SW types i s based on methods similar to ones d escrib ed in many papers (see 
reviews by IWimmer-Schweingruber et al.l (120061 ) and iTsurutani et al.l (120061 ) 
and references therein) and basically agrees with results of other authors but 
in contrast with other similar studies, we used general set of level criteria for 
all SW types and made identification for each 1-hour point. 



During 1976-2009 there were 798 magnetic storms with Dst < —50 nT. There 
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were data gaps in several parameters of OMNI database for 334 (42%) from 
them (these storms are denoted as "Ind eterminate" (IND]) and interplane- 
tary drivers were found for 464 storms ( Yermolaev et al. . 2009bl ). Magnetic 
storm is considered to be connected with specific SW type if its onset is ob- 
served in 2 hours after beginning and during this SW type (2 hours interval 
is the average d elay between appear a nce o f southward IMF and rea ction of 
magnetosphere ( Gonzalez and Echer . 2005 : Yermolaev et al.l . 2007al bl). THe 
statistics of magnetic storm distribution over different types of SW is fol- 
lowing: 145 storms were induced by CIR, 62 - MC, 16 1 - Ejec ta, 96 -Sheath 



(12 before MC and 84 before Eiecta) (lYermolaev et al.l . l2009bl ). In our previ- 



ous papers (lYermolaev et al.l . l2006l . l2007al JblJd) we considered magnetic storms 
with Dst < —60 nT and statistics of storms was a bit different (total number 
of storms is 623). 



As has been indicated above the d uration of main phase of magnetic storms 



varie s in a wide range o f 2-15 hours ([Vichare et al.l . l2005l : iGonzalez and Echer 



20051 : lYermolaev et al.l. l2006l. l2007bl'). ayerage durati on in the interval 1976- 



2000 is 7 ± 4 hours (lYermolaev et al.l . l2006l . l2007br ). We used SEA with 2 



reference time instants and located all onsets at "0" time of epoch and all Dst 
minima at "6" time. The times before onset {t < "0" time) and after Dst min- 
imum {t > "6" time) are real, but the time between onset and Dst minimum 
was re-scaled (proportionally increased/decreased). After this transformation 
all storms have equal durations of main phase in the new time reference frame. 
This interval of main phase has 5 sub-intervals. The time in the " 0" - " 6" in- 
terval for 2/3 storms was changed not more than by 1/3 of its duration. This 
"double" SEA method allows us to simultaneously study interplanetary con- 
ditions resulting in beginning and end of magnetic storms as well as dynamics 
(temporal variations) of parameters for storms with different durations. 



Taking into account that typical durations of large-scale SW types are sig- 
nificant ly shorter than the f ull duration of m a gnetic storms (see , for ex- 



ample. ICrooker et all (119991 ): iGoshng and Pzzol ([19991); iLvnch et all (I2003D : 
Lepping et al.l tOO^ ) and in the average they are 9 ± 4, 28 ± 12 and 20 ± 8 



hours for Sheath, ICME and CIR, respectively (lYermolaev et al.l . l2007al ). we 



restricted durations of curves calculated using double SEA by time from -12 
up to +24 hours. It is important that statistics decreases from main phase 
of storms to edges of (-12, +24) interval (especially for Sheaths) and errors 
may increase at the edges of interval (Standard deviations for different pa- 
rameters and different in terplanetary drivers are similar to data presented by 
(lYermolaev et al.l . l2007al )). In cases discussed below differences between the 



curves in figures are mathematically significant although sometimes they are 
less than corresponding standard deviations. In some another cases it is nec- 
essary to consider these differences as a tendency (hypothesis) rather than a 
proven physical fact. The further investigations are required to reduce this 
uncertainty. 
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3 Results 



Figures 1 and 2 present time variation of several interplanetary and mag- 
netospheric parameters for 798 magnetic storms with Dst < -50 nT during 
1976-2000, which were obtained by the double SEA method with 2 reference 
epoch zero times: Dst storm onset and Dst minimum (dashed lines which cross 
"0" and "6" in the time axis of figures, respectively) and for 6 interplanetary 
drivers: (1) MC, (2) Ejecta, (3) sum of MC and Ejecta, (4) Sheath (sum of 
Sheaths before MC and Ejecta), (5) CIR, and (6) IND ("Indeterminate") type 
of SW. Curves for different types of solar wind are presented by different sym- 
bols/color. Under figures designations and number of events for each SW type 
are specified, the number of points in a separate bin of curves in figures can 
be less than the specified number of events, especially at the interval edges. 
Figure 1 shows: (Left column) V - velocity, T - proton temperature, T/Texp 
- ratio of measured proton temperature to calculated t emperature Texp us- 



ing a verage dependence of temperature on velocity V (iLopez and Freeman 



19861 ). Dst index, (Right) n - density, Pt - thermal pressure, Pdyn - dynamic 



pressure, Kp index. Figure 2 presents other parameters: (Left column) Ey - 
electric field, Bz - GSM southward co mponents of IMF, Dst* - density cor- 
rected Dst index fiBurton et all . Il975h . AE index, (Right) By,Bx,B - GSM 
components and magnitude of IMF, j3 - ratio of proton thermal to magnetic 
pressure. 



For all interplanetary drivers of magnetic storms the onset starts in 1-2 hours 
after southward turning of IMF {Bz < 0) and the main phase stops in 1-2 
after disappearance of the southward IMF component. In new time scale in 
region of "0" - "6" there is no significant difference between behaviour of Dst 
(Dst*) index for different drivers. Nevertheless, it is possible to indicate a slight 
tendency that most sharp decrease of indices is observed for Sheath and MC 
while the largest values of southward component of IMF and electric field are 
in MC. The highest values of Kp and AE indices are generated by Sheath and 
MC. Slope of Dst (Dst*) index and values of Kp and AE indices for Ejecta and 
MC+Ejecta are less than for MC, and this fact is one of the most important 
reasons why we did not consider all ICME together and made selection of two 
subtypes of ICME: MC and Ejecta. The highest value of velocity V is observed 
in Sheath (difference relative to another curves is 70-100 km/s), temperature 
T and T/Texp (2-3 times) in Sheath and CIR, density n (1.5-2.0 times) in 
Sheath and CIR, thermal pressure Pt (5-7 times) in Sheath, dynamic pressure 
Pdyn (1.5-2.0 times) in Sheath and CIR, {3 (1.5-2.0 times) in Sheath and CIR, 
IMF magnitude B (2-5 nT) in Sheath, MC and CIR. There is no systematic 
difference in Bx and By components of IMF for different SW types in the 
region of "0" - "6" time. 

To study behaviour of interplanetary parameters separately during moderate 
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and strong storms we divide data presented in Figs.l and 2 in two groups 
with -100 < Dst < -50 nT (see Figs. 3 and 4) and Dst < -100 nT (Figs. 5 
and 6). This selection decreased data statistics in the new figures, especially 
for Figs. 5 and 6 (see number of events under figures) when number of points 
in a separate bin of curves at interval edges can be less than 10. Taking into 
account accuracy of our estimations, it is possible to tell that all specified above 
tendencies and features of relative behaviour of SW and IMF parameters and 
magnetospheric indices for different interplanetary drivers are the same for 
both storm sizes. 

One of important problems of connection between interplanetary conditions 
and magnetospheric processes is dependence of magnetospheric activity on 
temporal evolution of SW and IMF parameters including Bz and Ey. We 
found a consistency between time evolution of cause {Bz and Ey) and time 
evolution of effect {Dst, Dst*, Kp and AE indices) for the time interval of 
"0" - "6" using 3 procedures (see Figs. 7-9): (1. Left column in figures where 
abscissa is designated as Ey and Bz) simultaneous measurements, for ex- 
ample, dependence of Dst{f).vs.Ey{f),i = 0, ...,6, (2. Second column, ab- 
scissa is Ey{t — 1) and Bz{t — 1)) 1-h displaced measurements, for exam- 
ple, Dst{f).vs.Ey{f~^),i — 0, ...,6, and (3. Right column, abscissa is Ey{Y) 
and Bz{Y,)) dependence of indices on integral value of sources, for example, 
Dst\vs.Ey{Y.y = Jq Ey{T)dT = Eo%^^ = 0,...,6;k = 0,...,i. It is impor- 
tant to remind that all storms have equal durations in interval of "0" -"6". 

The left column of Fig. 7 shows that there is no any monotonic relation between 

interplanetary conditions and magnetospheric indices. If we take into account 
the 1-h delay between sources and effects (the second column), dependencies 
become more monotonic, i.e. there is delay between cause and effect. The right 
column of figure demonstrates that discussed processes have a " memory" and 
all dependencies are monotonic and almost linear. It is important to note that 
the most effective process of Kp and AE indices generation acts during Sheath, 
CIR and IND, the least effective process of Dst and Dst* generation - during 
MC. Additional selection of Fig. 7 data on storm size with -100 < Dst < -50 
nT and Dst < -100 nT (see Figs. 8 and 9, respectively) shows that main 
properties specified in Fig. 7 remain in the new figures. 



4 Discussion and Conclusions 

We performed an analysis of interplanetary conditions for 798 magnetic storms 
with Dst < -50 nT for the period 1976-2000 on the basis of the OMNI archive 
data. Our analysis has two special features. (1) Taking into account importance 
of epoch time selection we used the method of the double superposed epoch 
analysis including simultaneously two reference times: the onset of storm and 
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the minimum Dst index. (2) Taking into account different reaction of magneto- 
sphere on various interplanetary disturbances we categorized large-scale types 
of solar wind as interplanetary drivers of storms: CIR, Sheath, MC, Ejecta, 
ICME (e.g. both MC and Ejecta) and "Indeterminate" type. This methodical 
approach showed the following. 

First of all we would like to note that our method reproduced some well- 
known results and, in particularly, showed that independently on types of 
interplanetary drivers, the onset of magnetic storms begins in 1-2 hours after 
southward turn of IMF and recovery phase of storms begins in 1-2 hours after 
disappearance of this component of IMF. This is a good verification of used 
method. At the same time the method allowed us to obtain some new results. 



Various types of interplanetary drivers of magnetic storms have significantly 
different parameters. Particularly, Sheath and CIR (with respect to ICME) 
have higher density, dynamic and thermal pressures, temperature, /3-parameter 
as well as higher variance of the same parameters and magnitude and com- 
ponents of IMF. These differences can be significant, for example, dynamic 
pressure is larger by a factor of 1.5 - 2. Our an alysis confirms high importance 

of Sheath in generation of magnetic storms (IHuttunen and Koskinenl. |2004J: 

Yermolaev et al. . 2006 ; Huttunen et al. . 2006 : Yermolaev et al. . 2007a b^: Pulkkinen et al. 



2007al ) and indicates specific geoeffective conditions in Sheath. It is important 



to note that there are serious differences in parameters of MC and Ejecta 
(including their geoeffectiveness) , while differences in Sheath before MC and 
Sheath before Ejecta are not significant. This fact should be taken into account 
during analysis. 



Using double SEA method (transformation of main phases of storms to equal 
duration) allowed us to compare storms with different durations of main phase. 
In the re-scaled temporal reference frame the Dst. vs. time (Dst*. vs. time) de- 
pendencies for storms induced by different types of interplanetary drivers are 
close to each other and have approximately linear shapes in time interval of " 0" 
- "6". Comparison of Dst. vs. time and Dst*. vs. time dependencies for storms 
induced by Sheath showed that dynamic pressure Pdyn results in parallel 
displacement of Dst (Dst*). vs. time dependencies but does not change slope of 
these dependencies. Though our tr ansformation should mask the shorter main 
phase for Sheath- induced storms ( Yermolaev et al. . 2007bl : Pulkkinen et al. 



2007al ). our results speak in favor that the Dst. vs. time (Dst*. vs. time) depen- 



dencies for storms induced by Sheath are more abrupt. The Kp. vs. time and 
AE. vs. time dependencies are nonlinear and higher for Sheath- than for CIR- 
and ICME-induced storms. It means that these indices are generated by other 
mechanisms than Dst and Dst* indices. 



Dependencies of Dst (or Dst*) on the integral of Bz over time are almost lin- 
ear and parallel for different types of drivers. This fact can be considered as 
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an indication that time evolution of main phase of storms depends not only 
on current values of Bz and Ey but also on their prehistory. The differences 
between these lines are relatively small {ADst < 20 nT), nevertheless there is 
tendency that the MC curve lies higher (i.e. at equal values of integral of Bz 
the storm is smaller). Dependencies on integral of Ey are similar. Dependen- 
cies of Kp (and AE) on integral of Bz (and Ey) over time are nonlinear and 
nonparallel. The differences between these lines are relatively small {AKp < 1 
and AAE < 50 nT), nevertheless there is tendency that the Sheath and CIR 
curves lie higher than for another drivers. 
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Table 1 

List of results on interplanetary conditions resulting in magnetic storms obtained 
by superposed epoch analysis 



N 


Number (Years) 


Zero time 




Selection 


SW and IMF 


Reference 


1 


538(1963-1991) 


onset 




No 


B,Bx,By,Bz,V,T,n,Pdyn 


Tavlor et al. f 1994) 


2 


120(1979-1984) 


min Dst 




No 


Bz, n, V 


Maltsev et al. (1996) 


3 


150(1963-1987) 


turning Bz 




No 


Bz,Pdyn 


Davis et al. ("1997) 


4 


305(1983-1991) 


onset 




No 


Bz, Pdyn 


Yokovama and Kamide Q997) 


5 


1085(1957-1993) 


min Dst 




Dst 


Bz, Pdyn 


Loewe and Prolss ('1997') 


6 


130(1966-2000) 


onset 




No 


B,Bx,ByBz,|Bx[, \By\, 
|B2|,V,n,Pdyn 


Lvatskv and Tan (2003) 


7 


623(1976-2000) 


onset and min 


Dst 


SW types" 


B,Bx,By,Bz,V,T,n,Pdyn, 
nkT, (3, T/Texp 


Yermolaev et al. (2005. 2006. 2007a.b) 


8 


78(1996-2004) 


min Dst 




SW types'' 


B,Bz,dB/B,V,T,n, 


Mivoshi and Kataoka (2005) 


9 


549(1974-2002) 


min Dst 




Yes= 


B,Bx,ByBz,|Bx|, \By\, 
|B2|,Bs,VBs,V,n,T,Pdyn 


Zhane et al. (2006) 


10 


623(1976-2000) 


onset 




SW types" 


aB, aV, aT, an 


Yermolaev et al. (2007c) 


11 


28(19972002) 


onset and min 


Dst 


SW types'* 


Bz, Pdyn, V, Ey 


Pulkkinen et al. (2007a) 


12 


10(2004) 


onset 




No 


Bx,ByBz,B,e,V,n,Pdyn, 
A'/A.Ey 


Pulkkinen et al. (2007b) 


13 


29(1999-2002) 


onset, main phase, 


No 


Bz,Pdyn 








min Dst 








Hie et al. (2008) 


a 


(1) cm, (2) Sheath and (3) MC; - 


(1) CIR and (2) MC (Sheath + MC body); 





- (1) moderate storm at solar minimum, (2) moderate storm at solar maximum, (3) strong storm 
at solar minimum, and (4) strong storm at solar maximum; - (1) MC and (2) Sheath. 
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Fig. 1. Behavior of interplanetary parameters and magnetospheric Dst and Kp in- 
dices for magnetic storms with Dst < -50 nT generated by different interplanetary 
drivers: (1) all ICME (MC+Ejecta), (2) MC, (3) Ejecta, (4) Sheath, (5) CIR, and 
(6) "Indeterminate" (see designations in bottom of figure) during 1976-2000 on the 
basis of OMNI database obtained by double superposed epoch analysis method 
with two reference times: onset ("0" time, 1st dashed line) and minimum Dst index 
("6" time, 2nd dashed line). Presented parameters: (Left column) V - velocity, T 
- proton temperature, T/Texp - ratio of measured proton temperature to calcu- 
lated temperature Texp u sing average dependence of temperature on velocity V 
( Lopez and Freemanl . 19861 ). Dst index, (Right) n - density, Pt - thermal pressure, 
Pdyn - dynamic pressure, Kp index. 
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ALL storms 




Fig. 2. The same as in Fig.l for parameters: (Left column) Ey - electric field, 
Bz - GSM southwa rd components of IMF, Dst* - density corrected Dst index 



([Burton et all Il975l ). AE index, (Right) By, Bx, B - GSM components and magni- 



tude of IMF, P - ratio of thermal to magnetic pressure. 
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storms with -100<Dst<-50 nT 




Fig. 3. The same as in Fig.l for moderate storms with -100 < Dst < -50 nT. 
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Fig. 4. The same as in Fig. 2 for moderate storms with -100 < Dst < -50 nT. 
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Fig. 5. The same as in Fig.l for strong storms with Dst < -100 nT. 
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storms with Dst<=-100 nT 




Fig. 6. The same as in Fig. 2 for strong storms with Dst < -100 nT. 
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Fig. 7. Time evolution of Dst, Dst*, Kp and AE indices during time evolution of Bz 
and Ey for time interval of "0" - "6" using 3 procedures (see text): (1. Left column) 
simultaneous measurements, (2. Second column) 1-h displaced measurements, and 
(3. Right column) dependence of indices on integral value of sources. 
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Fig. 8. The same as in Fig. 7 for moderate storms with -100 < Dst < -50 nT. 
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Fig. 9. The same as in Fig. 7 for strong storms with Dst < -100 nT. 
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